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ABSTRACT 



We use the redshift-magnitude relation, as derived by Dabrowski 
(1995), for the two exact non-uniform pressure spherically symmetric 
Stephani universes with the observer positioned at the center of 
symmetry, to test the agreement of these models with recent observations 
of high redshift type la supernovae (SNIa), as reported in Perlmutter 
et al. (1997). By a particular choice of model parameters, we show 
that these models can give an excellent fit to the observed redshifts and 
(corrected) B band apparent magnitudes of the Perlmutter et al. data, 
but for an age of the Universe which is typically about two Gyr - and may 
be more than three Gyr - greater than in the corresponding Friedmann 
model, for which non-negative values of the deceleration parameter 
appear to be favoured by the data. We show that this age increase is 
obtained for a wide range of the non-uniform pressure parameters of the 
Stephani models. We claim this paper is the first attempt to compare 
inhomogeneous models of the universe with real astronomical data. 

Several recent calibrations of the Hubble parameter, from the Hubble 
diagram of SNIa and other distance indicators, indicate a value of 
H ~ 65, and a Hubble time of ~ 15 Gyr. Based on this value for H Q 
and assuming A > 0, the P97 data would imply a Friedmann age of at 
most 13 Gyr and in fact a best-fit (for q = 0.5) age of only 10 Gyr. Our 
Stephani models, on the other hand, can give a good fit to the P97 data 
with an age of up to 15 Gyr. The Stephani models considered here could, 
therefore, significantly alleviate the conflict between recent cosmological 
and astrophysical age predictions. The choice of model parameters is 
quite robust: in order to obtain a good fit to the current P97 data, 
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one requires only that the non-uniform pressure parameter, a, in one of 
the models is negative and satisfies \a\ < 3 km 2 s~ 2 Mpc -1 . This limit 
gives a value for the acceleration scalar, u, of order \u\ <,0.66 x lCT 10 r 
Mpc -1 , where r is the radial coordinate in the model. Thus, although 
the pressure is not zero at the center of symmetry, r = 0, the effect of 
acceleration is non-detectable at the center since the acceleration scalar 
vanishes there. However, the effect of the non-uniform pressure on the 
redshift-magnitude relation is clearly seen since neighbouring galaxies are 
not situated at the center and they necessarily experience acceleration. 
By allowing slightly larger, negative, values of a one may 'fine tune' the 
model to give an even better fit to the P97 data. 

Subject headings: Cosmology - age of the Universe - supernovae - 
relativity 
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1. Introduction 

The standard isotropic Friedmann cosmological models have naturally been 
the most widely investigated models in studies of the large-scale structure of the 
Universe. This is hardly surprising, in view of their mathematical simplicity and 
their generic prediction of an approximately linear Hubble expansion at low redshift, 
which is in excellent agreement with observational data (of. Strauss & Willick 
1995; Postman 1997). Even in Friedmann models, however, the relation between 
apparent magnitude and log redshift is in general non-linear at higher redshift and 
depends explicitly on the spatial curvature of the Universe - or equivalently on the 
deceleration parameter, q . 

For several decades astronomers have attempted to use the Hubble diagram 
of some suitable 'standard candle' (e.g. first-ranked cluster galaxies) to place 
constraints on the global geometry of the universe by comparing the observed 
redshift-magnitude relation of the standard candle with that predicted in Friedmann 
models with different values of go ( c -f- Peach 1970; Gunn & Oke 1975; Schneider, 
Gunn & Hoessel 1983; Sandage 1988). Results from such analyses have thus far 
proved inconclusive, however. Due to the intrinsic dispersion in the luminosity 
function of the standard candles available, one previously had to reach at least 
z ~ 1 before the predictions of models with different values of q became sufficiently 
distinct to be detectable; at the same time, however, the effects of luminosity and 
number density evolution also become important at these redshifts, and are very 
difficult to correct for. The situation for the Hubble diagram of quasars is equally - 
if not even more - problematical. Tonry (1993) suggested that the contraints on q 
from such studies were no better than — 1 < q < 1. 

Recently, however, it has been suggested that type la supernovae (hereafter 
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SNIa) represent a standard candle of sufficiently small dispersion to allow meaningful 
estimates of g now to be derived from the SNIa Hubble diagram at more moderate 
redshift. In Perlmutter et al. (1997; hereafter P97) a preliminary analysis is 
presented of seven distant SNIa in the range 0.35 < z < 0.50. A comparison of the 
SNIa magnitudes and redshifts with the predicted relation for various Friedmann 
models appears to exclude large negative values of go, and is best fitted by values 
close to go = 0.5. This poses a potentially serious problem for Friedmann models. 
Since many recent determinations of the Hubble constant (including a number of 
analyses using SNIa) suggest that H lies in the range 65 - 70, this would imply an 
age of the universe of less than 10 Gyr in the 'standard' f2 = 1,A = scenario. 
This result would appear to be in sharp conflict with recent astrophysical age 
determinations from e.g. globular clusters and white dwarf cooling (c.f. Chaboyer 
1995; Hendry & Tayler 1996) - a conflict which is only slightly alleviated by revisions 
to globular cluster age estimates in the light of results from the HIPPARCOS satellite 
(Chaboyer et al. 1997). Reducing the value of f2 lessens the conflict somewhat, but 
agreement is still only marginal if one accepts a robust lower bound for the matter 
density of f2 m = 0.3, as has been suggested by several different methods of analysing 
large-scale galaxy redshift surveys (c.f. Strauss & Willick 1995). This situation has 
helped to give a renewed impetus to models with a positive cosmological constant 
(c.f. Liddle et al. 1996) which contributes an additional component, Q\, to make up 
the critical density and at the same time extends the age of the Universe by up to 2 
Gyr - depending on the value of H and fl m . However, because of the relation 

go = T^m - ^A, (1) 

it is clear that a positive value of g is incompatible with a positive value of A unless 
the matter density is at least two-thirds of the critical density. The go = case, 
assuming fl m = 2/3, = 1/3 and H = 65, would give an age of the Universe of 
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just over 11 Gyr; as go increases the age is decreased still further. Thus, if the results 
of P97 prove to be correct and the deceleration parameter is non-negative, then 
the conflict between cosmological and astrophysical age predictions remains firmly 
unresolved - at least if H >, 65. Independent results showing that a positive value 
of A is incompatible with the so-called VLBI data (Kellermann 1993), using the 
angular diameter test, were obtained by Krauss & Schramm (1993) and Stelmach 
(1994). 

In this paper we propose one method to alleviate this age conflict by considering 
some inhomogeneous cosmological models in which the relation between the age 
of the Universe and a generalised Hubble constant is more general than in the 
Friedmann case. Despite some theoretical plots of the observational quantities for 
inhomogeneous models (e.g. Goicoechea & Martin-Mirones 1987, Moffat & Tatarski 
1995, Dabrowski 1995, Humphreys, Maartens & Matravers 1997) this paper is - as 
far as we are aware - the first to compare these models with real astronomical data. 
In particular, we show that taking an inhomogeneous model into account allows us 
to obtain a good fit between the predicted redshift-magnitude relation and P97 data, 
but for an age of the Universe which is several Gyr older than in the Friedmann 
case. The models under consideration have been discussed before and are known 
as Stephani Universes (c.f. Kramer et al. 1980; Krasihski 1983; Dabrowski 1993). 
In these models the energy density g depends just on the cosmic time, similarly to 
the Friedmann models, but the pressure, p, is a function of both spatial coordinates 
and a time coordinate; hence the models are usually referred to as 'inhomogeneous 
pressure universes'. In the spherically symmetric case under consideration, the 
pressure is just a function of time and radial coordinate which means that its 
values are the same on spheres, r = constant, around the center of symmmetry but 
differ from sphere to sphere. This essentially means that there is a spatial pressure 
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gradient and particles are accelerated in the direction from high-pressure regions to 
low-pressure regions. This effect is usually described by the acceleration vector ii r 
which in the case of spherical symmetry has only one (radial) component, or by the 
acceleration scalar u (cf. Eqs. 2.18 and 2.23 of Dabrowski 1995). The acceleration 
represents the combined effect of gravitational and inertial forces on the fluid which, 
in fact, similarly as in Newtonian physics, are unable to be separated. As a first 
approximation we assume the observer is placed at the center of symmetry, which 
results in no pressure gradient at the observer's postion. This, in a sense, contradicts 
the Copernican Principle, but can easily be overcome by applying the formulae for 
a non-centrally observer given in Section 5 of Dabrowski (1995) - an appropriate 
generalisation once larger samples of high-redshift supernovae become available. 

In the standard approach we neglect the effect of pressure (i.e. we take 
pressureless dust - with pressure p = 0) as we evaluate chaotic velocities of galaxies 
to be small. This results in taking the acceleration, ii, also to be zero (cf. Ellis, 1971, 
for a discussion of the relation between these quantities). However, if there was large 
flux of neutrinos or gravitational waves for instance, this assumption would not be 
correct and we would need to take radiation pressure (p — ^g) into account. This has 
been of course investigated for isotropic cosmologies (Dabrowski & Stelmach 1986, 
1987) and all the observational quantities have been found. Our main point here 
is, however, that early universe processes such as phase transitions (for details see 
Vilenkin 1985, Kolb & Turner 1990, Vilenkin & Shellard 1994) may result in having 
different exotic types of matter (e.g. cosmic strings) with many different types of 
equations of state. In the easiest case (straight cosmic strings) they may end up 
with the exotic equation of state p = —\q 1 but in general, the equation of state can 
be more complicated (e.g. Vilenkin & Shellard 1994 and de Vega & Sanchez 1994, 
in the context of superstrings) or spatially dependent (e.g. Narlikar, Pecker & Vigier 



- 8- 



1991). The latter would especially be the case of our interest. Of course the standard 
energy conditions of Hawking and Penrose might be violated (cf. Hawking & Ellis 
1973) which also happens for inflationary models for instance and our considerations 
here are, in a sense, on the same footing as those phenomena. As for the Stephani 
models, which do not admit any global barotropic equation of state, it has been 
shown that there exists a consistent nonbarotropic equation of state and the full 
thermodynamical scheme exists (Quevedo & Sussman 1995, Krasihski, Quevedo & 
Sussman 1997). 

Regardless of the physical background of the models under consideration, one 
of our main tasks here is to draw attention to the entire class of inhomogeneous 
models which could be a useful alternative to Friedmann models in helping to 
resolve the apparent incompatibility of measurements of the Friedmann cosmological 
parameters. Even if the final outcome (after a thorough comparison with data) 
shows that the universe is indeed isotropic and homogeneous this conclusion must 
be drawn by applying some 'averaging scale-dependent procedures' ( cf. Ellis 1984, 
Buchert 1997) since we evidently cannot see the universe to be like that on smaller 
scales. Being spherically symmetric, Stephani models can also be applied to a local 
underdense/overdense spherical region embedded in a globally isotropic Friedmann 
universe (Moffat & Tatarski 1995) in some analogy to the so-called 'Swiss Cheese' 
model (Kantowski, Vaughan & Branch 1995). 

The reader interested in more generic models should be referred to the recent 
review by Krasihski (1997), as well as to some earlier papers concerning the most 
popular generalization of the Friedmann models such as the spherically symmetric 
Tolman Universes which are inhomogeneous density pressure-free dust shells 
(Tolman 1934; Bondi 1947; Bonnor 1974). Their properties have been studied 
quite thoroughly in Hellaby & Lake (1984, 1985) and Hellaby (1987, 1988) and 
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the observational relations for Tolman models were studied by Goicoechea & 
Martin-Mirones (1987), Moffat & Tatarski (1995) and quite recently by Humphreys, 
Maartens & Matravers (1997). However, in none of these cases has a comparison 
with real astronomical data been carried out. 

The outline of this paper is as follows. In Section 2 we reproduce the 
redshift-magnitude relations for the two Stephani models considered here, as recently 
derived in Dabrowski (1995). In Section 3 we briefly describe the SNIa data of P97. 
In Section 4 we fit these data to the redshift-magnitude relations of both Friedmann 
and Stephani models and thus obtain best-fit values for the model parameters. We 
then discuss the results of these fits and compare the age of the Universe given 
by the best-fit model parameters in the Friedmann and Stephani cases. Finally in 
Section 5 we summarise our conclusions. 

2. The Redshift-magnitude relation for inhomogeneous pressure models 

Recently Dabrowski (1995) has considered the redshift-magnitude relation for 
Stephani universes. Two exact cases were presented and the predicted relations were 
plotted for a range of different parameter values. The relations were defined following 
the method of Kristian & Sachs (1966), of expanding all relativistic quantities in 
power series and truncating at a suitable order. Approximate formulae, to first order 
in redshift z, for Model I and Model II respectively were given by 




+ 1.086 l + 4a 



(arg + 6r + d) 
(2ar + b) 2 



z 



(2) 



and 



m B = M B + 25 - 5 log 10 h 5 log 10 cz 



o 
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+ 1.086 ^ + ^0^)^ (3) 

which are essentially equations (5.6) and (5.10) respectively of Dabrowski (1995). 
[] Here m-Q and Mb denote apparent and absolute magnitude respectively in the 
B band, r denotes the current age of the Universe and c = 3 x 10 5 km s _1 is the 
velocity of light. The constants a, h and d are parameters of Model I and a is a 
parameter of Model II. Convenient units for these parameters are: [a] = km 2 s -2 
Mpc -1 , [b] = km s -1 , [d] = Mpc and [ac 2 ] = (km s" 1 Mpc)~t. From the definition 
of the acceleration scalar (c.f. Dabrowski 1995) we conclude that the parameters 
which relate directly to the non-uniformity of the pressure are a in Model I and a in 
Model II. In Model I, b plays a similar role to the coefficient of time in the expression 
for the scale factor (R oc r v where p is any power) in Friedmann models and can be 
considered as an immanently Friedmannian parameter of the Stephani models, while 
a and a are completely non-Friedmannian. 

The models considered here are spherically symmetric, which means that we 
can have both centrally placed and non-centrally placed observers. For simplicity 
the redshift-magnitude relations reproduced above correspond to a centrally placed 
observer. Dabrowski (1995) also derived relations for the case of a non-centrally 
placed observer which, although more general, introduced several additional free 
parameters. The main difference in this more general case is that the apparent 
magnitude depends on the position of the source in the sky, and renders comparison 

x But written now for B band observations, which are the case under study in this 
paper. Note also that we have corrected a typographical error which appeared in 
Dabrowski (1995) - the factor of 4 in the final term of Eq. (2) replacing the factor of 
2 in Eq. (5.6) of that paper. We thank Chris Clarkson for drawing our attention to 
this error. 
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with the Friedmann case more complicated. We thus consider only centrally placed 
observers in this paper. 

Note that these formulae are truncated at first order in z and thus would become 
increasingly inaccurate if applied to redshifts greater than or equal to unity. Since 
the redshifts of the P97 preliminary data extend only to z ~ 0.5 we proceed with 
the first-order expressions, and will also use the first-order expression for Friedmann 
models when comparing the fits. In future work, as the database of high redshift 
SNIa grows, we will extend the approximate redshift-magnitude relations to higher 
order, as required. 

The reader is referred to Dabrowski (1995) for a detailed discussion of the 
derivation of the above formulae. Note, however, that for Model I the expression for 
the generalised scale factor, R(t), as a function of cosmic time, r, is given by (Eq. 
2.11 of Dabrowski 1995) 



However, R(r) does not have to be positive (Dabrowski 1993, Section IV.A and Fig. 
6) for the Stephani models. For the subclass under consideration, R(t) easily relates 
to the spatial curvature of the models 



and the curvature index is not constant in time as for Friedmann models. In 
principle, one can restrict R(t) to be positive (this is especially reasonable, if we 
want to obtain the full Friedmann limit), which ends up with the simple relation for 
the spatial curvature of the models being positive for negative non-uniform pressure 
parameter a and negative for positive a. Since R(t) = at the singularity (the 
Big-Bang) we can require that R(t) — > as r — > (i.e. we set the origin of our time 
coordinate at the Big Bang) and thus demand that d is identically zero. Therefore, 



R(t) = ar 2 + br + d. 



(4) 



k{r) = -A- 2 R(r), 



(5) 
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according to condition (2.13) of Dabrowski (1995), which, in fact, allows one to have 
the Friedmann limit for the Stephani models under consideration, we have 

b 2 = 1, (6) 

for the values of b. Without loss of generality we assume that b = +1 (cf. the 
discussion above about the meaning of b in Friedmann models) - leaving only one 
free parameter of the Model I, which is a. 

Other important physical quantities of Model I are as follows: the energy 
density 

8nG . . 3 

C( r ) = T^TTT^> (7) 



the pressure 



c 4 KV y (aT 2 + br) 2 
8ttG , , 1 



and the acceleration scalar 



u = -1%r. (9) 



c 2 



From the above one can see that the finite density singularities of pressure appear at 
r — > oo where there is the antipodal center of symmetry. We assume that we are 
placed at the center of symmetry at r = so we have these singularities far away 
from us (see Fig. 6 of Dabrowski 1993). Of course we cannot live at the singularity 
of pressure. 



2 Note that in Stephani models we use the so-called isotropic radial coordinate 
which results in the Friedmann metric to be taken in the isotropic form. It relates to 
the commonly used nonisotropic coordinate via the relation = r/(l + (l/4)r 2 ) 
(for details see Section II of Dabrowski 1993 ). Because of isotropy at every point one 
cannot differentiate the centers r = and r = oo in a Friedmann universe, which is 
not the case for the Stephani models. 
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At the center of symmetry the fluid fulfils the barotropic equation of state 
p = —^g (as for straight cosmic strings, cf. Vilenkin 1985), while at r — > oo the 
pressure goes to either plus or minus infinity. Then, assuming R(t) > 0, it diverges 
to minus infinity if a > 0, and to plus infinity if a < 0. In such a case, the particles 
are accelerated away from a high pressure region at r = to low pressure regions 
at r 7^ 0, if a > 0, and toward a low pressure region at r = from high pressure 
regions at r ^ 0, if a < 0. Of course if R(t) < the situation is the opposite. The 
acceleration scalar is zero at r = and it diverges at r — > oo. 

In the case of Model II the time-dependent curvature index is given by (8 plays 
here the same 'Friedmannian' role as b in Model II) 

fc( r ) = -a/3 2 c - M, (10) 

while the energy density, pressure and acceleration scalar are given respectively by 
(Dabrowski 1993, Appendix C) 



8ttG _ I 1 3n 

2 

T3 



C 4 £00 = jj^ --a > (11) 



87rG 2a 4 a/5 2 . , ,, , 

— o-p(t) = — - o ^ r + a r ' ( 12 ) 

" 2 73 3 T 3 



and 



■it = —-a/3r. (13) 



3. The SNIa observations of Perlmutter et al. (1997) 

Type la supernovae are thought to be the result of the thermonuclear disruption 
of a white dwarf star which has accreted sufficient matter from a binary companion to 
reach the Chandrasekhar mass limit. For several decades they have been considered 
as suitable (nearly) standard candles for the testing of cosmological models because 
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of the relatively small dispersion of their luminosity function at maximum light and 
the fact that they are observable at very great distances. In recent years the Hubble 
diagram of SNIa has been used by a number of authors to obtain estimates of the 
Hubble constant (c.f. Riess, Press & Kirshner 1996; Hamuy et al. 1995, 1996; Branch 
et al. 1996) and the motion of the Local Group (Riess, Press & Kirshner 1995). 
P97 consider the redshift-magnitude relation of SNIa at high redshift, observed by 
the 'Supernova cosmology project', as a means of constraining q . In P97 SNIa are 
not treated as precise standard candles, but a 'stretch factor' correction is applied 
to account for the correlation between SNIa luminosity and the shape of their light 
curve. 

In this paper we use the redshifts and B band magnitudes - with and without 
light curve shape corrections - as presented in Table 1 of P97, to which the reader 
is referred for details of their observing strategy, data reduction procedures and 
magnitude error estimates. 

4. Comparison of the data with Friedmann and Stephani Models: 

results and discussion 

4.1. Friedmann Models 

Figure 4 of P97 shows the Hubble diagram of their SNIa compared with 
the theoretical magnitude-redshift relations for a Friedmann model with different 
combinations of Q m and fl\. While P97 argue correctly that one should generally 
express the Friedmann magnitude-redshift relations in terms of f2 m and 
separately, and not just in terms of their combination via q , for the redshift range 
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of the P97 data one may adequately approximate the relation by 



m B = Mb + 5 log 10 cz + 1.086(1 - q )z, 



(14) 



where 



Mb = M B - 5 log 10 H + 25, 



(15) 



with the corresponding expression for the corrected B band magnitudes. For reasons 
which will become clear when we consider the Stephani models, it is useful for us 
to write Eq. (14) in this form, in terms of go- Since we will make use of similar 
expressions for the Stephani universes, we construct for the Friedmann case the 
(reduced) chi-squared statistic 



x 2 = 



n 



i i=i . 



m 



pred , 



a(i) 



(16) 



where n is the number of SNIa, m B hs (i) and a(i) are respectively the observed B 
band apparent magnitude and error estimate of the i th SNIa, and m B rcd (i) is the 
predicted B band apparent magnitude of the i th SNIa, for a given value of q , derived 
from equation 5 (or its equivalent for the corrected magnitudes). Following P97 we 
adopt Mb = -3.17 ± 0.03 and A^ B ,corr = -3.32 ± 0.05. 

From Eqs. (14) and (16) it follows that q , the maximum likelihood (equivalently 
minimum y 2 ) estimate of qo, is given by:- 



E 



where 



and 



Xi = 1.086z(i) 



(17) 



(18) 



Vi = m B bs (i) - Mb -5\ogcz(i) - 1.086z(i). 



(19) 
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Substituting the appropriate values from P97 we find that q = 0.48 for the 
uncorrected magnitudes and % = 0.50 for the corrected magnitudes. Thus we see 
that applying the magnitude corrections has negligible effect on the best-fit value of 
go for the P97 data. 

Figure 1 shows the value of x 2 \ f° r —0.5 < go < 1, both for the uncorrected 
(solid line) and corrected data (dashed line). We can see from this Figure that the 
corrected data consistently give a slightly better fit to a Friedmann model than 
the uncorrected data for a wide range of values of q . Both the uncorrected and 
corrected data give an acceptable fit over the range < go < 1 (and in fact over 
a somewhat wider range for the corrected data) but both give a poor fit for large 
negative values of go- 

Table 1 quantifies the goodness of fit of the SNIa data to a number of Friedmann 
models with different values of the cosmological parameters fl m , and the age of 
the universe, to - all for a Hubble constant H = 65 km s _1 Mpc -1 . Column (3) 
shows the corresponding value of go, calculated from Eq. (1). Column (5) gives 
the reduced \ 2 °f t ne fit to all seven SNIa, while column (6) gives the reduced x 2 
obtained using the five SNIa with corrected magnitudes. 

EDITOR: PLACE TABLE |Tj HERE. 

It is clear from Table 1 that one cannot obtain, with Ho ~ 65 km s" 1 Mpc -1 , 
an acceptable fit to either the corrected or uncorrected data and at the same time 
ensure an age of the universe in excess of 14 Gyr. We discuss the situation for other 
values of H below. 
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4.2. Stephani Model II 



We now compare the SNIa data with the predicted magnitude-redshift relations 
of the Stephani models. We consider first Model II, and the relation given by Eq. 
(3). If we compare Eqs. (3) and (14) we see that in the limit as z — > these 
equations are identical if and only if 



In other words for nearby SNIa, Stephani Model II predicts the same linear 
redshift-magnitude relation as do Friedmann models, and with an age of the universe 
equal to two thirds times the inverse of the Friedmann Hubble constant. This is 
precisely the age of a Friedmann universe which is flat with a zero cosmological 
constant. In particular, if H ~ 65 km s _1 Mpc -1 then independent of the value 
of the parameter, a, the age of the universe r in Model II is approximately 10 
Gyr, which certainly appears to be too low to be consistent with astrophysical age 
determinations. Hence it would seem that Model II is not particularly useful in 
resolving the current age conflict since the age is inextricably linked to the value of 
the Friedmann Hubble constant: as soon as the latter is specified then so too is the 
age of Model II. 

The link between the magnitude-redshift relation for Model II and the 
Friedmann case is, nonetheless, interesting for the following reason. Note that Eq. 
(3) may be rewritten as 



(20) 



m B = 



M B + 25 - 5 log 10 — + 5 log 10 cz 



o 



+ 



1.086(1 -q )z, 



(21) 



where 




(22) 
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which means that for any given age of the universe, To, we can choose the parameter, 
a, so that the magnitude-redshift relation for Model II is identical in form to Eq. 
(14), with t = ^Hq 1 . The crucial difference is that, whereas in the Friedmann case 
with A = 0, Eq. (20) implies that q = 0.5, in the Stephani case we still retain the 
freedom to specify a relation which is equivalent to any value of g by suitable choice 
of a. 

In particular, by choosing a < one can obtain a magnitude-redshift relation 
which corresponds to a Friedmann model with g > 0.5. This is in full analogy 
to Friedmann models if the relation for the curvature of the Stephani models is 
taken into account (Dabrowski 1995, Eq. 2.14). It shows that the time dependent 
curvature index (Eq. 10) for a < is positive (if cosmic time, r > 0), while for a > 
it is negative. The pressure (Eq. 12) is positive or negative for a being positive or 
negative respectively at the center of symmetry r = and it diverges to either plus or 
minus infinity (depending on the values of other parameters) at the antipodal center 
of symmetry r — > oo and the particles are either accelerated away or towards r = 0. 
Bearing in mind the effect of curvature (Eq. 10) of the models one can roughly say 
that the inclusion of non-uniform pressure mimics a flat Friedman model go — 1/2 
to become curved - positively curved for a < 0, and negatively curved for a > 0. 

While in the Friedmann case g > 0.5 would imply an age of the universe 
t < ^Hq 1 , in the case of Model II we still have r = §# _1 . Model II would, 
therefore, be of considerable interest if SNIa (or other) observations were to suggest 
that g > 0.5, which certainly cannot be ruled out on the basis of the P97 data alone. 
As an illustrative (if somewhat extreme) example, consider the case where fl m = 2 
and Q\ = 0, so that g = 1 for the Friedmann model. As can be seen from Eq. (22) 
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the Friedmann and Model II magnitude-redshift relations are identical when 



Whereas the age of the Friedmann model with q = 1 would be reduced by about 



unlikely, in view of a variety of other observations of large scale structure and CMBR 
anisotropies, this serves as an interesting example of how the Stephani models can be 
compatible with high redshift observations over a larger region of parameter space 
than Friedmann models. 



One of the reasons why Model II is not particularly useful as an extension of 
the Friedmann case is that the effect of the non-uniform pressure (manifest via the 
parameter a) only becomes apparent at high redshift. The situation with Stephani 
Model I is different, however. We can see from Eq. (2) that the effect on the 
magnitude-redshift relation of the non-uniform pressure parameter a is immediate. 
In particular, therefore, even at low z Model I does not in general reduce trivially to 
a specific Friedmann case. 

Note that after setting the parameter b from Model I equal to unity (provided 
d — 0) we can rewrite Eq. (2) to depend only upon the non-uniform pressure 
parameter a. Thus 
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15% compared with the Einstein de Sitter age (i.e. To — 0.57 H 1 ), for the Stephani 
Model II we still have r = ^Hq 1 . Although the scenario of go > 0.5 appears highly 



4.3. Stephani Model I 



m B = M B + 25 + 5 log 10 cz 




+ 1.086 1 +4 



(ar 2 + r ) 



z. 



(24) 



(2ar + l) 2 



As a means of estimating what range of values of a and tq will give an acceptable 
fit to the P97 data it is useful to note further that we may recast Eq. (24) in the 
form 



m B = M B + 25 + 5 log 10 cz - 5 log 10 H 



+ 1.086(1 - q )z, 



(25) 



where 



H = 



2ar + 1 



(26) 



and 



q = -4a 



ar 2 + T Q 



(27) 



(2ar + l) 2 ' 



Eq. (24) now takes the same functional form as Eq. (14), as was similarly 
pointed out in Dabrowski (1995), with H and % replacing H and q . We can 
think of Hq (which is one third of the expansion scalar 6 of the model) and q as a 
generalised Hubble parameter and deceleration parameter which are related to the 
age of the universe in a different way from the Friedmann case. The key question of 
interest here is therefore whether one can construct generalised parameters, H and 
q , which are in good agreement with the P97 data but which correspond to a value 
of t which exceeds that Friedmann age with H = H and q = %. The fact that we 
can write the Model I redshift-magnitude relation in the form of Eq. (25) confirms, 
however, that our choices of r and a are certainly not arbitrary. Combinations 
of r and a which give a large negative value of qo, for example, would clearly be 
incompatible with the SNIa Hubble diagram - just as was the case for Friedmann 
models with q « 0. 

In order to estimate the parameters r and a we construct the reduced 
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chi-squared statistic:- 

(2« 



/ = — ^- y 



t=l L 



where n is the number of SNIa and mg red (i; r , a) is obtained for the i th SNIa from 
Eqs. (25), (26) and (27). We determine Mb from Hamuy et al. (1996), adopting 
their best-fit value of Hq = 63.1 km s _1 Mpc -1 determined from four local calibrating 
SNIa. Thus, 

M B = .Mb + 5 log # -25 = -19.17 ±0.03 (29) 



and 



M BiCorr = A^Bxorr + 5 log H - 25 = -19.32 ± 0.05. (30) 



The dependence of Eq. (28) on r and a is non-linear, making a plot of the 
surface z = x 2 ( T o, a ) difficult to interpret. We therefore consider slices through this 
surface. Moreover for plots of x 2 at constant tq it is useful to plot \ 2 as a function of 
a -1 . Figure 2 shows x 2 as a function of a" 1 for r = 13 Gyr, using the uncorrected 
P97 data. The behaviour of x 2 is seen to be rapidly varying for values of a -1 around 
zero, but is essentially flat for all |a _1 | >,0.3. Thus, provided |a| < 3 km 2 s~ 2 Mpc^ 1 , 
we see that the goodness of fit of Model I to the data is essentially independent of 
the value of a, and depends only on tq. A very similar curve is obtained for the 
corrected magnitudes. 

We can understand the rapidly varying behaviour of x 2 f° r small values of | cz 1 1 
by considering the behaviour of H and go in Eqs. (26) and (27). We see that when 
a -1 = — 2r we have Hq = and \qo\ — > oo 0, so that x 2 ~> 00 • It therefore follows 
that for r = 13 Gyr = 13.26 x 10~ 3 s Mpc km -1 a singular value of x 2 occurs when 



3 This is a special situation, and moreover one which contradicts the observations 
since it would mean that the present day is exactly the turning point of the cosmic 
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a -1 ~ —0.026, and x 2 varies very rapidly close to this value. However, the range 
of very small a^ 1 is not of interest to us since it deviates too far from Friedmann 
models. 

Figures 3a and 3b show plots of x 2 as a function of a" 1 but now with t = 15 
Gyr, for the corrected and uncorrected magnitudes respectively. A narrower range 
of values of a' 1 is shown, in order to better illustrate the behaviour of x 2 for small 
| a 1 1 . We find that x 2 again tends to infinity when a -1 = — 2tq, and is again 
essentially flat for all |a _1 | >, 0.3. Note also that the asymptotic value of x 2 is a little 
smaller than that for r = 13 Gyr, and moreover that there exists a narrow range of 
values of a for which x 2 dips appreciably below its asymptotic value. 

Figures 4a and 4b, on the other hand, show plots of x 2 as a function of r for 
a = —1.0, for the uncorrected and corrected magnitudes respectively. This value of 
a is chosen to be representative of the asymptotic behaviour of x 2 \ essentially the 
same plots would be obtained for all \a\ <, 3 km 2 s~ 2 Mpc -1 . We can see that Model 
I gives a good fit to the data for r in the range 13 to 15 Gyr. 

Figures 5a and 5b show the values of H and % respectively as a function of 
t , again for the representative value of a = —1.0. Also shown for comparison are 
the best-fit values of H = 63.1 km s" 1 Mpc -1 and q = 0.5, obtained from Hamuy 
et al. (1996) and section 4.1 respectively. We see from Figure 5b that % is almost 
independent of r over the range shown, increasing from q ~ 0.04 (r = 10 Gyr) to 
g — 0.08 (to = 20 Gyr). The dependence of H on r is rather more pronounced, 

evolution. Following Eqs. (7) and (8) the energy density at the moment t = — l/2a 
is d>-nG/c A Q = 48a 2 and the pressure 8nGp = g(—l/3 + l/6r 2 ). In the Friedman limit 
a^Owe have r — > oo, and p = g = which is an everlasting, empty, flat and static 
universe - a subcase of little or no physical interest. 
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however: good agreement with the Hamuy et al. value is found in the age range 14 
to 16 Gyr. 

The behaviour of H and % in Figures 5a and 5b makes sense when we consider 
the form of Eqs. (26) and (27) for ar << 1. To first order in ar these reduce to 

H = -(l + ar ) (31) 

To 

and 

g = -4ar . (32) 
Thus we see that as ar — > 0, H — > Tq 1 and g — > 0. 

The potential usefulness of Model I is now apparent. In the limit where ar — > 0, 
the age of the universe in this model is increased by 50% compared with the Einstein 
de Sitter age, giving for example r = 15 Gyr (compared with only 10 Gyr) for 
Hq ~ 65. Of course when ar — > 0, g — > also, so that the more meaningful 
comparison is that between Model I and a Friedmann model with q = 0. Taking 
fi m = 0.3, go = and H = 65 gives a Friedmann age of only 12.5 Gyr, however, so 
that the Model I age is still 2.5 Gyr greater than the Friedmann age. Taking larger 
values for the matter density results in a bigger difference between the Friedmann 
and Model I ages. 

Figures 5a and 5b are indicative of the limiting behaviour of H and go for 
small ar . As we already remarked for the x 2 plots, one obtains similar plots for all 
\a\ <, 3 km 2 s~ 2 Mpc -1 , with H strongly varying as a function of r but g much 
more weakly dependent on r . The smaller the value of \a\ the closer g lies to zero 
at given r - as is obvious from Eq. (32). The sign of a does have some bearing on 
the goodness of fit of the model to the P97 data, however. Although it can be seen 
from Figure 1 and Table 1 that the current data give an acceptable fit for g = 0, the 
fit rapidly deteriorates for negative values of go. Thus, if a > (i.e. a high pressure 
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region at r = away from which particles are accelerated) then go —> from below, 
and a value of a = 3 km 2 s~ 2 Mpc -1 would imply g ~ —0.2 for r = 15 Gyr, which 
gives only a marginally acceptable fit to the P97 data. If a < 0, on the other hand 
(i.e. a low pressure region at r = towards which particles are accelerated) then one 
can obtain much better fits: e.g. a = — 3 km 2 s~ 2 Mpc -1 implies go — 0.2 for r = 15 
Gyr. 

We have emphasised the limiting behaviour of Model I, for small a, in order to 
make clear that the usefulness of the model is a fairly robust result and is not too 
sensitive to the exact value of a which is chosen - although it is true that negative 
values of a are favoured. It is particularly noteworthy that the limit \a\ <3, may 
be considered as the restriction on this parameter from the observational data and 
allows us not to be too far from the range where Friedmann models are valid. In 
other words one can obtain a good fit to the P97 data, with a significantly larger 
age, but without requiring that Model I departs too much from a Friedmann model. 

According to Eq. (9), \a\ <3 km 2 s~ 2 Mpc— 1 translates to a limit on the value 
of acceleration scalar of <,0.66 x 10~ 10 r Mpc -1 , where r is the radial coordinate 
in the model. Notice that, although the pressure is different from zero at the center 
r = the effect of acceleration is not detectable at the center since the acceleration 
scalar (and of course vector) vanishes there. 

If we allow slightly larger negative values for the non-uniform pressure 
parameter, a, then by careful choice of a and r we can obtain fits which give 
significant positive values of g , while retaining a significant difference between 
the Stephani and Friedmann ages. Although these fits require slightly more 'fine 
tuning' they are clearly in much closer agreement with the best-fit Friedmann value 
of go = 0.5. Some examples of fits of this type are given in Table 2. The final two 



-25 - 



columns of Table 2 show the age, rp, of the universe in a Friedmann model with 
A = and with Q m + Q A = 1 respectively. 

EDITOR: PLACE TABLE g HERE. 

Some general trends are evident from Table 2. Note that in all cases we see that 
as To increases and \a\ becomes smaller, the values of Hq and go ar e both reduced 
and the goodness of fit to the corrected and uncorrected data gradually deteriorates. 
For t > 16 Gyr, the goodness of fit quickly becomes unacceptably large: although 
by suitable choice of a one can ensure that go remains positive, the value of H also 
reduces and overall the fit deteriorates. Further decrease in \a\, for r > 16 Gyr, 
increases the value of H , but pushes g closer to zero, so that the goodness of fit 
remains poor. This behaviour can also be easily seen from Eqs. (31) and (32). 

It would seem, therefore, that an age of To = 15 — 16 Gyr represents the 
upper age limit from Model I with the P97 data - at least if one adopts the SNIa 
calibration of H . Moreover, if subsequent analysis of larger samples of SNIa serve to 
tighten the limits on a positive value of %, then this limiting age could perhaps be 
reduced to r ~ 14 Gyr. The important point to note, however, is that in this case 
the age limits on Friedmann models would also be reduced. As can be seen from 
Table 2, a value of go ~ 0.5 can be well fitted by Model I with tq = 14 Gyr, which 
still represents an increase in the age of the universe of more than 3 Gyr compared 
with the go = 0.5 Friedmann model with either zero cosmological constant or critical 
density. 
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5. Conclusions 

In this paper we have considered the two exact non-uniform pressure spherically 
symmetric Stephani universes, and have compared the redshift-magnitude relations 
derived for these models in Dabrowski (1995) with the recent SNIa observations of 
P97. We have investigated the extent to which, by suitable choice of the Stephani 
model parameters, we may obtain good fits of the P97 data to the predicted 
redshift-magnitude relations but for universes which are older than their Friedmann 
counterparts. 

We emphasize that we have considered only the case of centrally placed 
observers, which results in having zero pressure gradient in our location. Although 
this is clearly a special case, it is mathematically the simplest possibility and merits 
consideration first. It can be extended relatively simply to the case of a non-centrally 
placed observer using the formulae given in Section 5 of Dabrowski (1995). However, 
such a generalisation introduces additional model parameters which make apparent 
magnitude a function of both redshift and direction the sky. In principle we could 
have estimated these parameters in this paper, but the small size of the P97 sample 
would make such a parameter fit statistically meaningless. Indeed, even for the case 
of a centrally placed observer, ideally one should consider a much larger supernovae 
sample. We will consider the more general case in future, when the number of 
observed supernovae has significantly increased. 

We have found that the age of the universe in Stephani Model II is, in fact, 
independent of the non- uniform pressure parameter, a, and is equal to the age of 
an Einstein de Sitter Friedmann model, i.e. r = ^Hq 1 . This model would be of 
considerable interest if the total density of the universe were greater than the critical 
density, since the age of the corresponding Friedmann model would then be less than 
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the Einstein de Sitter age. Since there exists no compelling observational evidence 
to suggest that the universe is closed, however, Model II is of limited use as it would 
in general predict an age which was smaller than its Friedmann counterpart with the 
same value of the Hubble constant. 

We have shown that Stephani Model I would be of considerably greater interest, 
however. We have found that the redshift-magnitude relation predicted for Model I 
can be expressed in terms of two parameters: the age, To, of the universe and the 
non-uniform pressure parameter a. One can write the redshift-magnitude relation 
in exactly the same form as in the Friedmann case, introducing an effective Hubble 
parameter, Hq, and deceleration parameter, qo, which are non-linear functions of a 
and r . We have shown that for a wide range of different values of a we can obtain 
good fits to the P97 data for a universe of age up to 15 Gyr, which is typically 
two or three Gyr greater than the corresponding Friedmann model. These fits are 
quite robust, requiring only that \a\ <, 3 km 2 s~ 2 Mpc -1 , which gives the value of 
the acceleration scalar u of the order <, 0.66 x lCT 10 r Mpc -1 , where r is the 
radial coordinate of the model. Then, although the pressure is not zero at the center 
of symmetry r = 0, the effect of acceleration is non-detectable at the center since 
the acceleration scalar vanishes there. However, this effect is easily extracted from 
the redshift-magnitude relation since neighbouring galaxies are not situated at the 
center and they necessarily experience acceleration. 

The above robust fits are for the limiting case where the product ar is small, 
and imply an effective deceleration parameter, q , close to zero - a value which is 
certainly not as yet ruled out by the P97 data, although the fit to a value of q~ = 0.5 
is currently somewhat better. By some fine-tuning of the Model I parameters, one 
can obtain good fits with q ~ 0.5 and r ~ 14 Gyr. While an age of only 14 
Gyr may still be in conflict with astrophysical age determinations, the conflict is 
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considerably worse for Friedmann models: the age of an Hq ~ 65, go ~ 0.5 critical 
density Friedmann universe is only ~ 10 Gyr, and for closed models with q ~ 0.5 
the age is even smaller. 

Thus, we find that Model I can give an age of the universe which is consistenly 
and robustly between two and three Gyr older than the oldest acceptable open or 
flat Friedmann models. 

Since the preliminary results of P97 were first presented, there have been 
several important developments in the measurement of fundamental cosmological 
parameters. The recalibration of the RR Lyrae distance scale has revised age 
estimates of the oldest globular clusters to t = 11.7 ± 1.4 Gyr (c.f. Chaboyer et 
al. 1997). This undoubtedly lessens the conflict with the standard (Q = 1, A = 0) 
cosmological model - particularly if one argues for a value of Hq ~ 55 (c.f. Tammann 
1996). If one requires a 'gestation period' of around 1 Gyr between the Big Bang 
and the formation of the first globular clusters, however, then agreement with the 
standard model is still only marginal - even for Hq = 55 - and open Friedmann 
models would appear to be favoured. Since we have argued in this paper that the 
P97 data does not yet exclude models with q ~ 0, it is only fair to point out 
that open Friedmann models with A ^ 0, q ~ and Hq ~ 55 offer a comfortable 
resolution of the age problem, in the light of the Chaboyer et al. results. Adopting, 
for example, Q m = 0.5, = 0.25 and H = 55 gives a Friedmann age of rp = 14.0 
Gyr. 

It is important to recognise, however, that this agreement rests crucially upon 
the value of Hq. If, instead, one adopts the most recent estimate of H from the 
HST distance scale Key Project: H = 73 ± 6 ± 8 km s" 1 Mpc" 1 (Freedman 1996), 
the above Friedmann age reduces to only 10.6 Gyr, and for the standard model 
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(with q = 0.5) is only 8.9 Gyr. Moreover, the impact on H of the HIPPARCOS 
recalibration of the LMC distance modulus has recently been shown by Madore 
and Freedman (1997) to be less significant than had previously been reported (c.f. 
Feast & Catchpole 1997). Thus it would seem that rumours of the end of the age 
problem are perhaps somewhat premature. If, indeed, the value of Hq lies close 
to that obtained by the HST Key Project then we note that the Stephani models 
considered here can still give an age of up to approximately 12.5 Gyr, with q ~ 0.5, 
and 13.4 Gyr, with q ~ 0. While the data clearly do not yet present a case for 
abandoning Friedmann models, equally they do not rule out the possible need to 
do so in the future - when bounds on H and go are tightened - and the Stephani 
models investigated here could indeed prove to be very important. 

In this paper we have considered only a particular class of inhomogeneous 
models, in order to illustrate their potential usefulness in addressing the apparent 
conflict between the observed values of the Friedmann model parameters. In future 
work we will extend our treatment to a wider class of non-Friedmann models 
and test their compatibility with the Hubble diagram of high-redshift objects and 
other cosmological observations. Such a comparison will be greatly enhanced by 
having larger samples of distant SNIa - a development which the modern observing 
strategies adopted by P97 and other groups will shortly provide. 

As for the further progress in our comparison of the Stephani models with 
astronomical data one should of course investigate such standard tests as galaxy 
number counts, the angular size-redshift relation or microwave background 
anisotropies as far as to adopt the second order terms in redshift z 2 . This can 
indeed be done after some relatively tedious calculations, again using the powerful 
power series methods originally given by Kristian & Sachs (1966, see also Ellis 1971 
for more detailed discussion). Many of these issues have been studied previously 
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for Friedmann models (e.g. Dabrowski & Stelmach 1986, 1987; Stelmach, Byrka & 
Dabrowski 1990) and will be the subject of future work. 
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Fig. 1. — The value of x 2 -, for —0.5 < q < 1, both for the uncorrected (solid line) 
and corrected data (dashed line). 

Fig. 2. — Plot of the variation of x 2 with a -1 , for a fixed value of r = 12 Gyr, obtained 
from a comparison of the apparent magnitudes predicted by Stephani Model I with 
the uncorrected magnitudes of P97. 

Fig. 3. — Plot of the variation of x 2 with a^ 1 , for a fixed value of r = 15 Gyr, obtained 
from a comparison of the apparent magnitudes predicted by Stephani Model I with 
the magnitudes of P97. Figure 3a is for the uncorrected magnitudes and Figure 3b is 
for the corrected magnitudes. 

Fig. 4. — Plot of the variation of x 2 with r , for a fixed value of a = —1.0, obtained 
from a comparison of the apparent magnitudes predicted by Stephani Model I with 
the magnitudes of P97. Figure 4a is for the uncorrected magnitudes and Figure 4b is 
for the corrected magnitudes. 

Fig. 5. — Plot of the effective Friedmann parameters, H and go, as a function of r , 
for a fixed value of a = —1.0. Dotted lines indicate the best-fit Friedmann values for 
H (Figure 5a) from Hamuy et al. (1996) and go (Figure 5b) from section 4.1. 
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Table 1. Goodness of fit, expressed as a reduced x 2 , of the SNIa data to 
Friedmann models with different combinations of Q m and Q\. Column (3) gives the 

corresponding value of qo and column (4) the age, to, of the universe assuming 
H = 65 km s" 1 Mpc -1 . Columns (5) and (6) respectively give the reduced x 2 of the 
fit to all seven SNIa and to the five SNIa with corrected magnitudes. 
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Table 2. Fits of Stephani Model I universes with significantly positive values of q 
to the P97 data. The final two columns show the age, r F , of the universe in a 
Friedmann model with the same values of H and q , with A = and with 

fl m + Q A = 1 respectively. 
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